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INVESTIGATION  OF  TRANSONIC  GAS  FLOWS 


O.M.  Belotserkovskll  * 


This  report  considers  some  problems  of  transonic  gasdynamlcs  related 
to  the  theory  of  planar  and  three-dimensional  flows  of  a perfect  gas  which 
have  been  worked  out  mainly  by  the  author  and  his  collaborators.  Results 
of  analytical  investigations  of  a number  of  exact  properties  of  solutions 
which  describe  the  flow  around  bodies  with  a detached  shock,  obtained  by 

E. G.  Shifrin,  O.M.  Belotserkovskii , et  al.,  are  presented  without  proof  in 
the  form  of  theorems  (proofs  are  to  be  found  In  the  papers  cited).  Results 
of  numerical  solutions  (due  to  O.M.  Belotserkovskii,  Yu.  M.  Davidov, 

F. D.  Popov,  et  al.),  obtained  by  the  methods  of  integral  relations  and 
"large  particles",  are  presented  for  the  complete  problem  of  transonic  flow 
around  blunt  bodies  whose  profiles  possess  sharp  corners.  The  boundaries 
of  distinct  minimal  domains  of  influence  in  mixed  flow  about  blunt  bodies 
are  also  examined. 


Professor,  Member  of  the  Academy  of  Sciences  of  the  USSR,  Computing 
Center,  Academy  of  Sciences,  40  Vavilova  St.,  Moscow  V-333,  U.S.S.R. 

ii  . 


PREFACE 


Selected  results  of  this  report  and  its  companion  "New  Computational 
Models  in  Continuum  Mechanics,"  were  communicated  in  seminar  lectures  given 
by  Prof.  Belotserkovskii  at  several  American  universities  during  a 4-week 
visit  in  November  - December,  1978.  In  response  to  numerous  requests  Prof. 
Belotserkovskii  made  available  the  manuscripts  of  two  reports  which  describe 
in  some  detail  the  computational  techniques  employed  in  the  numerical  solu- 
tions of  the  problems  surveyed  In  his  seminar  talks. 

This  report  surveys  results  of  research  at  the  Computational  Center  of 
the  USSR  Academy  of  Sciences,  Moscow,  over  the  past  decade  in  some  problems 
of  planar  and  three-dimensional  transonic  flows  of  a perfect  gas.  Although 
most  of  these  results  have  previously  appeared  in  various  Soviet  journals, 
this  report  nevertheless  gives  a coherent  review  of  the  advances  in  the 
computation  of  transonic  flows  at  one  of  the  foremost  centers  of  the  Soviet 
Union. 

The  hallmarks  of  their  numerical  techniques  are  that  they  are: 

(1)  typically  differentially  and  globally  conservative,  and 

(2)  careful  attention  is  made  to  minimize  truncation  errors  while  at 
the  same  time  the  favorable  properties  of  the  computational  schemes 
are  exploited,  e.g.,  the  effective  viscosity  of  the  finite-difference 
equations  (to  promote  calculatlonal  stability)  and  (physically 
interpretted)  the  ability  to  compute  the  essential  features  of  se- 
parated regions  of  recirculating  flows  (e.g.,  wakes)  wholey  within 
the  framework  of  the  Euler  equations. 

This  work  Involves  a careful  and  mutually  beneficial  blending  of  analysis 
and  numerical  computations  to  develop  a complete  picture  of  complex  flows. 


iii 


The  numerical  methods  are  Illustrated  by  a great  variety  of  computational 
results  which  are  supported  by  analysis  and  compared  with  experimental  re- 
sults In  some  cases.  Wherever  possible  English  translations  of  the  refer- 
ences have  been  cited. 

The  careful  typing  of  the  edited  manuscript  by  Miss  Vicki  Brewer  de- 
serves a special  note  of  appreciation.  Finally,  I am  pleased  to  acknowledge 
Contract N00014-79-M-0022  from  the  U.S.  Office  of  Naval  Research  which  made 
possible  the  publication  and  distribution  of  this  report. 

W.  L.  Melnik,  Editor 

Professor,  Aerospace  Engineering  Dept. 

University  of  Maryland  at  College  Park 
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1.  INTRODUCTION 


Much  attention  has  recently  been  given  to  the  investigation  of  transonl 
problems  of  gasdynamics.  Such  flows  occur,  for  example,  in  nozzles,  around 
a body  with  a detached  bow  shock,  around  convex  corners,  as  well  as  on  wings 
with  imbedded  local  supersonic  regions. 

Difficulties  arise  in  theoretical  analysis  of  such  problems  primarily 
because  of  the  fact  that  there  is  no  complete  mathematical  formulation  of 
the  problem  in  a number  of  cases.  In  mixed  flow  domains  complex  phenomena 
occur:  so-called  secondary  "suspended"  compression  shocks  can  originate  be- 
hind the  bow  wave  and  for  supercritical  flight  Mach  number  local  supersonic 
domains  can  also  be  closed  off  by  compression  shocks.  The  mechanism  by 
which  such  phenomena  occur  has  not  been  sufficiently  studied  and  questions 
regarding  the  region  of  influence  and  its  adaptation  to  a change  in  flow 
parameters  also  remain  open.  Classical  analytical  methods,  developed  for 
plane  potential  flows,  which  predominantly  utilize  the  hodograph  plane,  do 
not  as  a rule  apply  here.  Numerical  solutions  of  these  problems  involve 
many  difficulties,  and  require  the  construction  of  special  schemes  for  the 
integration  of  elliptic-hyperbolic  equations. 

If  there  is  a break  (corner)  in  the  generator  of  a blunt  body  within 
the  domain  of  Influence  of  the  flow,  and  the  free  stream  is  altered  so  that 
the  local  speed  of  sound  occurs  at  this  corner,  then  the  analysis  of  such  a 
flow  is  complicated  substantially  by  the  existence  of  this  singularity.  The 
flow  expansion  around  the  corner  will  be  in  the  mixed  transonic  regime  and 
Is  characterized  by  an  abrupt  change  in  both  magnitude  and  direction  of  velo 
city.  Moreover,  a secondary  compression  shock  can  occur  In  the  supersonic 
zone  of  this  flow  which  substantially  affects  the  whole  flow  patten.  It 
should  also  be  noted  that  the  computational  results  for  the  flow  in  the 


domain  of  influence  of  the  bluntness  (M-domain)  serve  as  the  initial  data 
for  the  calculation  of  the  supersonic  zone.  Construction  of  the  solutions 
in  the  transonic  domain  for  a body  with  a corner  must  be  made  with  special 
care  since  even  slight  inaccuracies  in  the  calculations  will  not  permit  con- 
tinuation of  the  computation  into  the  supersonic  domain. 

A study  of  the  flow  properties  around  blunt  bodies  at  low  supersonic 
Mach  number  is  of  no  less  interest.  As  the  free  stream  Mach  number  decreases, 
the  domain  of  mixed  flow  influence  increases.  In  examining  such  a problem 
it  is  necessary  to  take  account  of  the  transonic  nature  of  the  flow  in  the 
zone  between  the  sonic  line  and  the  limiting  characteristic  separating  the 
M-domain.  Perturbations  in  the  transonic  domain  affect  the  shape  of  the 
sonic  line,  and,  therefore,  the  whole  flow  in  the  mixed  zone.  The  solution 
of  the  boundary  value  problem  becomes  increasingly  sensative  to  changes  in 
initial  data,  so  that  round-off  errors  increase,  and  an  instability  in  the 
computation  is  manifested.  All  this  demands  the  construction  of  special 
numerical  schemes.  Consequently  it  is  important  to  identify  the  existence 
of  boundaries  of  distinct  minimal  domains  of  influence  of  the  bluntness. 

Only  careful  numerical  experiments  utilizing  high-speed  computers  give 
quantitative  data  and  a complete  picture  of  the  flow  in  such  complex  problems. 
An  attempt  is  made  to  analyze  and  discuss  some  transonic  flow  properties 
(the  solutions  in  the  neighborhood  of  the  sonic  line,  the  formation  of  se- 
condary compression  shocks,  the  shape  of  the  minimal  domain  of  influence, 
etc.)  from  the  aspect  of  analytical  and  numerical  solutions. 

It  should  be  mentioned  that  precisely  by  numerical  means  are  secondary 
compression  shocks,  non-monotone  sonic  lines,  etc.,  successfully  detected 
and  constructed.  At  the  same  time,  the  conditions  for  spoiling  the  continu- 
ous solution,  and  the  regularities  associated  with  different  kinds  of  domains 
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of  Influence,  have  been  delineated  analytically.  In  examining  the  flow  j 

around  a corner,  the  effective  construction  of  numerical  schemes  turned 

out  to  be  possible  only  when  utilizing  asymptotic  methods  of  solution  of 

the  differential  equations.  Therefore,  the  combination  of  both  analysis 

and  numerical  solutions  prove  to  be  quite  fruitful. 

A number  of  papers  have  been  devoted  to  the  study  of  these  questions.  | 

We  consider  in  this  report  the  results  of  investigations  of  the  direct 

two-dimensional  (planar  or  axisymmetric)  problem  obtained  by  using  the  > 

methods  of  integral  realtions  [17, 23, 24, 26, 28, 29, 36]  and  "large  particles" 

[31,32,37],  which  represents  part  of  an  ongoing  research  program  in  tran- 
sonic flows. 
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2.  HODOGRAPH  PLANE 


L 


In  analyzing  plane  stationary  flows  behind  a shock  It  Is  convenient  to 
use  the  gasdynamics  equations  written  In  a local  streamline  coordinate  sys- 


tem: 


O . m2)  - m2  g- , Mae  . -k«  #■ 
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or  In  the  form  [1] 
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Here  X is  the  velocity  magnitude,  made  dimensionless  by  the  maximum 
adiabatic  velocity,  M - the  Mach  number,  B - the  slope  of  the  velocity  vec- 
tor measured  counter-clockwise,  p - the  pressure,  pQ  - the  total  pressure, 
and  k is  the  adiabatic  index;  3/3S-J , 3/3s2  - derivatives  with  respect  to 
the  directions  of  the  velocity  vector  and  its  normal  obtained  by  rotation 
of  the  velocity  vector  counter-clockwise  through  an  angle  ir/2. 


Equations  (1)  and  (2)  can  be  transformed  (by  using  differential  geo- 
metry formulas)  into  equations  in  a coordinate  system,  where  y ,<f>  are 
families  of  streamlines  and  their  orthogonal  trajectories  since  h^3/3s^  = 

= 3/3$ , h23/3s2  = q 3/34' ; the  Lame  coefficients  h1  ,h2  are  given  by  the  for- 
mulas [7] 


l dAnPn 

hjXexp[-  / — 2^-]  = *(<»),  h2q(x)  = F(f); 

M 

,(0-a<J$L-  ^ x2),/(k*1>  . 


(3) 


Here  $(x),  Fff)  are  arbitrary  functions  governing  the  numbering  of  the  lines 
in  the  families. 

In  the  case  of  irrotational  flows  (for  pQ  * const),  the  Chaplygin  equa- 
tions can  be  obtained  from  the  equations  in  coordinates  (by  a change  in 
the  role  of  the  dependent  and  independent  variables). 
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Equations  (3)  permit  the  Investigation  of  singularities  of  a mapping 
in  the  $,Y'  plane.  Indeed 


hence,  for  0 < A < (k+1  )ly^(k-l  )”1//2  a mapping  In  the  plane  Is  locally  one- 

-to-one  (x,y  are  Cartelsan  coordinates  In  the  physical  plane). 

Let  us  consider  the  mapping  of  the  velocity  (x,B)  and  pressure  (p,@)  In 
the  hodograph  plane;  the  former  Is  more  convenient  for  potential  flow  In- 
vestigations, and  the  latter  for  rotational  flow  investigations. 

It  follows  from  the  expressions  for  th^  Jacoblans  J = 3(p,B)/3(x,y) , 

I = 3(x,e)/3(x,y)  of  the  transformations,  using  (1)  and  (2),  that  the  mapping 
of  subsonic  domains  In  the  pB  plane  (and  in  the  AB  plane  in  the  potential 
flow  case)  can  only  have  Isolated  singularities  [2,3];  however,  this  does 
not  mean  that  the  Image  of  the  subsonic  domain  In  the  pB  plane  cannot  be 
many  sheeted  (on  the  whole);  the  flow  can  only  be  single-valued  In  the  physi- 
cal plane. 

A mapping  of  supersonic  domains  in  the  pB  plane  can  have  folds;  their 
edges  (on  which  J changes  sign)  are  called  branch  lines.  A branch  line  can 
be  a characteristic  If  J undergoes  a discontinuity  upon  passing  through  It; 

In  the  general  case  a branch  line  Is  the  envelope  of  characteristics  of  one 
family  In  the  pB  plane,  and  the  geometric  locus  of  cusps  of  the  second 
family  of  characteristics.  A mapping  of  a rotational  flow  In  the  XB  plane  can 
have  branch  lines  even  Irr  the  subsonic  domain.  In  the  general  case,  a branch 
Una  of  this  mapping  Is  a characteristic  only  In  tn-a  potential  flow  case  [3]. 
A classic  Illustration  of  the  flow  having  a fold  in  the  XB  plane  Is  a Laval 
nozzle;  the  branch  lines  are  here  characteristics  Issuing  downstream  from 
the  center-line  of  the  nozzle. 
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The  edge  of  the  fold  of  the  Inverse  mapping  Is  called  a limit  line.  The 
presence  of  a limit  line  Indicates  that  the  flow  whose  Image  has  been  found 
in  the  hodograph  plane  can  not  be  realized  in  the  physical  plane. 

The  characteristics  equation  in  the  pB  plane  have  the  form: 


± dB 


1,2 


(4) 


Hence  It  follows  that  for  1 < M < «°  the  first  family  characteristic  always 
makes  a negative  angle  with  the  p-axis,  and  a characteristic  of  the  second 
family  a positive  angle;  In  some  cases  this  permits  establishment  of  the 
form  of  the  branch  line  (If  It  exists). 

The  characteristics  In  the  XB  plane  are  given  by 

± dB,  o = - /(M2-l ) d«,nx  - d*nprt.  (4.) 


which  describe  epicycloids  in  the  case  of  potential  flow  (d«,npQ  * 0). 

THEOREM  1 . The  Image  In  the  XB  plane  of  a segment  of  characteristic 
contained  in  the  domain  of  monotonicity  of  the  function  pQ(v)  will  Intersect 
an  epicycloid  of  the  same  family  not  more  than  once.  A segment  of  a charac- 
teristic of  the  first  (second)  family,  drawn  In  the  direction  of  diminishing 
total  pressure  from  the  point  of  Intersection  with  the  epicycloid  of  the 
same  family.  Is  located  above  (below)  this  epicycloid  [4]. 

THEOREM  2.  The  Image  in  the  XB  plane  of  a characteristic  drawn  In 
the  domain  behind  the  shock  located  between  two  epicycloids  of  the  same 
family  which  are  separated  by  a distance  (in  the  B direction)  not  exceeding 

pom1n 


where  Pomax  and  Pom1n  are  the  maximum  and  minimum  values  of  the  total  pres' 
sure  behind  the  shock  (depending  on  the  free  stream  velocity)  and  k Is  the 
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adiabatic  Index. 


THEOREM  3.  If  the  velocity  vector  is  continuous  at  some  point  in  the 
physical  plane,  then  the  characteristic  at  the  corresponding  point  in  the  pe 
plane  Is  either  a smooth  curve,  or  has  a cusp  (the  latter  can  be  realized 
only  In  the  case  when  Intersection  with  a branch  line  occurs  at  this  point). 
This  property  of  the  characteristics  In  the  xp  plane  only  exist  under  the 
additional  condition  of  continuous  differentiability  of  the  total  pressure 
and  the  velocity  vector  at  the  corresponding  point  of  the  physical  plane. 
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3.  MINIMUM  DOMAIN  OF  INFLUENCE  OF  MIXED  FLOW 


The  minimum  domain  of  Influence  of  mixed  sub-  and  supersonic  flow  (the 
M-domaln)  Is  characterized  by  the  fact  that  small  perturbations  of  Its 
boundaries  are  propagated  throughout  this  domain.  In  the  general  case  the 
M-domain  boundary  is  delineated  by  segments  of  the  body-profile  which  may 
Include  a tangential  discontinuity  (corner),  the  axis  of  symmetry,  the 
bow  shock,  the  limiting  characteristic,  and  the  sonic  line. 


3.1  GENERAL  RESULTS 

In  order  to  establish  the  location  of  the  M-domaln  In  the  hodograph 
plane.  It  Is  necessary  to  know  the  sign  of  the  Jacobian  of  the  mapping  at 
points  of  the  boundary. 

THEOREM  4.  The  following  Inequalities  hold: 

1.  On  a straight-line  wall 

I * 3(x,B)/a(x,y)>  0,  I * 3(x,B)/3(x,y)  * 0 for  x < 1 ; 

2.  On  a free  boundary  J > 0 

3.  On  the  shock  polar  J > 0,  I $ 0 for  M < A(M-k), 

where  M^  Is  the  free-stream  Mach  number,  and  A Is  some  constant;  A | 1 for 
< M0U)t  where  MQ(k)  Is  some  constant. 

THEOREM  5.  The  Image  of  a convex  corner  In  the  pB  and  XB  planes  Is 
mapped  by  a characteristic  of  the  potential  flow  which  has  no  turning 
points*. 

THEOREM  6.  In  potential  flow,  the  Image  In  the  XB  plane  of  the  neigh 


* This  Is  proved  In  [5]  for  potential  flow. 
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borhood  of  a convex  corner  covered  by  characteristics  emerging  from  this 
point  Is  located  on  one  side  of  the  epicycloid  mapping  the  corner. 

Let  K denote  a point  of  the  sonic  line  at  which  It  Is  convex  relative 
to  the  lines  <}>  = constant  (sonic  line  at  the  point  K+  turns  the  convexity 
towards  the  supersonic  velocity  region,  and  at  K_  towards  the  subsonic 
domain).  Let  L denote  a point  of  the  sonic  line  at  which  the  sign  of  the 
streamline  curvature  changes  during  passage  through  It.  If  the  vortlcity 
Is  zero  at  the  point  K,  then  It  coincides  with  L. 

THEOREM  7.  There  are  no  points  K_  located  within  the  potential  flow 
domain. 

THEOREM  8.  If  there  Is  displacement  along  the  sonic  line  In  potential 
flow  so  that  the  subsonic  domain  lies  on  the  left,  the  velocity  vector  will 
be  rotated  monotonely  clockwise  [6]. 

THEOREM  9.  The  slope  of  the  velocity  vector  3 changes  monotonely  with 
displacement  along  a sonic  line  segment  In  a rotational  flow  which  does  not 
contain  the  points  K or  L.  A change  In  the  direction  of  variation  of  3 
occurs  at  each  of  the  points  K or  L [7]. 

Thus,  branch  lines  of  the  mapping  In  the  A3  plane  (an  odd  number  at 
each  of  the  points)  pass  through  the  points  K and  L. 

THEOREM  10.  If  the  acceleration  of  the  flow  at  the  sonic  line  Is 
bounded,  then  a(p,3)/  a(x,y)  - 0 at  the  point  L. 

Let  > denote  the  angle  through  which  It  Is  necessary  to  turn  the  tan- 
gent to  the  sonic  line  counter-clockwise  (at  the  sonic  point  of  the  shock) 
so  that  It  would  coincide  with  the  velocity  vector. 

THEOREM  1 1 . If  the  curvature  of  the  shock  wave  at  the  sonic  point 
does  not  vanish  or  become  Infinite,  then  y > ***  for  < N^(k)  and 
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Y < hit  for  > MQ( k)  [8]. 

The  constant  MQ(k)  Is  here  the  same  as  in  THEOREM  4 (see  Fig.  1 and  2). 

Theorem  11  can  be  generalized  partially  to  the  case  of  ax i symmetric 
flow.  Let  o denote  the  slope  of  the  shock  to  the  free  stream  velocity 
vector  measured  counter-clockwise. 

THEOREM  12.  If  the  curvature  of  the  shockwave  at  the  sonic  point  of 
an  axisymmetric  flow  does  not  become  zero  or  Infinite,  then  the  following 
inequalities  hold: 

1.  y < htt  for  a < 0 < and  M > M (k). 

00  0 

2.  y > hit  for  -hit  < a < 0 and  M^  < MQ(k). 

In  the  remaining  cases  the  angle  y (Fig.  2)  depends  on  the  shock  curvature, 
the  distance  from  the  axis  of  symmetry,  and  the  free  stream  velocity  [9], 

Let  6 denote  the  angle  between  the  sonic  line  and  the  profile  on  the 
downstream  side  of  the  sonic  line. 

THEOREM  13.  If  the  vorticity  is  zero  along  the  contour  of  a smooth 

profile,  then  0 < & < hit  on  the  convex  profile,  and  hit  < 6 < it  on  the  con- 

cave profile  (if  the  acceleration  8A/8s1  does  not  vanish  or  become  infi- 
nite at  the  sonic  point  of  the  profile)  [8]. 

The  next  theorem  follows  from  11  and  13. 

THEOREM  14.  If  the  assumptions  of  Theorems  11  and  13  are  satisfied. 

In  the  case  of  symmetric  flow  around  a smooth  convex  profile  with  a detached 
shock  an  odd  number  of  points  K exists  on  the  sonic  line  for  M^  > MQ(k) 

(at  least  one  point  K+);  for  Mm  < MQ(k)  the  number  of  points  K Is  zero  or 

even  (the  number  of  points  K+  equals  the  number  of  points  K_)  [8], 
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THEOREM  15.  At  a convex  corner  of  the  profile,  the  sonic  line  Is 
orthogonal  to  the  "subsonic"  direction  of  the  angle  (In  transonic  flow). 

This  property  has  been  proved  in  [10]  for  potential  flow. 

The  theorems  considered  above  assist  In  classifying  the  minimum  domains 
of  influence  for  an  unbounded  flow  around  a profile  (or  axlsymmetric  body) 
and  are  in  good  agreement  with  results  of  calculations. 

THEOREM  16.  In  symmetric  flow  around  a smooth  convex  profile  (If  the 
assumptions  of  Theorems  11  and  13  are  satisfied),  the  number  of  branch  lines 
of  a mapping  In  the  xe  plane  which  Intersect  the  sonic  line  Is  zero  or  is 
even  If  Mw  < MQ ( k ) ; the  number  of  these  branch  lines  Is  odd  If  Moo>  MQ(k)  [7] 

The  following  theorems  characterize  the  flow  In  the  neighborhood  of 
the  point  K of  the  sonic  line. 

THEOREM  17.  If  a point  K exists  on  the  sonic  line,  at  which  the  flow 
acceleration  does  not  become  zero  or  Infinite,  then  the  following  Inequali- 
ties are  satisfied  [7]: 

/ 9X  \Z  > 1 d£s  _ 1 /ds\2 

8sl  < Rklt+Ty  ^2  R2k2(k+1)  d<|> 

at  the  point  K+  or  K_.  Here  R Is  the  gas  constant,  and  s the  entropy. 

Let  AK+B  denote  the  line  formed  by  the  characteristics  Issuing  down- 
stream from  the  point  K+  (we  consider  the  direction  "downstream"  to  be  along 
the  streamlines  from  the  subsonic  domains). 

THEOREM  18.  At  least  one  branch  line  of  the  mapping  In  the  xe  plane 
which  passes  through  the  point  K+  In  the  supersonic  domain  will  be  down- 
stream of  the  line  AK+B  [7]. 


11 


If  the  flow  acceleration  at  the  point  K does  not  become  zero  cr 
infinite,  the  flow  in  the  neighborhood  of  this  point  can  be  investigated 
analytically  [11].  Equations  of  transonic  rotational  flow  are  derived  by 


the  small  parameter  method 


u -as-  - -2L  - 0 , -a-  + JBf  • «° , 

3$°  3<|>0  3$°  34<0 

“ s (k+i)(wrw02)  , W0  = , w,  =4^4.  . 

K 

B = ^ = s,2*\  - e3  , ♦ = £V,  r - £*°  . 


Here  e Is  a small  parameter  characterizing  the  size  of  the  neighborhood 
of  the  point  K. 

The  equations  obtained  have  an  exact  solution,  analogous  to  the 
exact  solution  of  the  transonic  potential  flow  equations,  which  was  used 
In  [12]  to  Investigate  the  flow  In  the  neighborhood  of  the  center  line 
of  a Laval  nozzle  under  the  condition  that  the  acceleration  does  not  be- 
come zero  or  Infinite. 

This  solution  has  the  form: 

2 2 
u * A4>  - *2,  V = (k+1  )W0<J>  + A2**  - Y3  , 


A « -(k+1)  £- 

3S1  K 


The  solution  so  obtained  can  be  analyzed.  For  example,  a three- 
-valued  mapping  function  can  be  constructed  In  the  hodograph  plane;  the 


geometry  of  the  branch  line  Is  studied  as  a function  of  the  parameter  o>, 
which  characterizes  the  velocity  (Fig.  3).  Also  the  characteristics 


2 

pattern  (Fig-  4)  Is  of  considerable  Interest,  particularly,  when  w > A 
(the  point  K_) . Analogous  Investigations  have  been  made  In  [13,14]  for 
axisymmetrlc  and  three-dimensional  rotational  flows. 


3,2  DOMAINS  OF  INFLUENCE  IN  FLOW  AROUND  BODIES 
WITH  A DETACHED  SHOCK  WAVE 


Let  us  examine  the  boundaries  of  various  kinds  of  minimal  domains  of 
Influence  In  the  flow  around  smooth  plane  or  axisymmetrlc  bodies  with  a 
detached  shock.  A schematic  diagram  of  the  domains  of  Influence  realized 
for  such  flows  in  given  In  Fig.  5 (the  velocity  vector  at  the  point  K is 
orthogonal  to  the  sonic  line). 

The  essential  characteristics  of  each  type  of  domain  of  influence  are 
the  slope  of  the  sonic  line  to  the  velocity  vector  on  the  shock  (y)  and  on 
the  body  (s).  The  minimum  domain  of  influence  shown  in  Fig.  5.  holds  for 
the  following  combinations  of  values  of  these  angles  [8,9,17]: 

1st  type  : y oblique  , 6 acute; 

2nd  type  : y acute  , <5  acute; 

3rd  type  : y acute  , 6 oblique. 

Since  the  angle  5 Is  always  acute  In  the  plane  case  (v  = 0)  In  the  flow 
around  a convex  profile  (Theorem  13),  the  third  type  of  domain  of  Influence 
Is  not  realized.  The  transition  from  the  first  into  the  second  type  occurs 
for  M * M„(k)  (Theorem  11).  The  graph  of  M (k)  Is  given  In  Fig.  2,  where 
MQ  = 1 .6358,1 .6895  and  1 .7421  for  k = 1.2, 1.4,  and  1.6  respectively.  Flow  patterns 
and  domains  of  Influence  for  M = 1,1.5,3,3.10  (k  = 1.4,  circular  cylinder), 

oo 

obtained  from  numerical  calculations,  are  presented  In  Fig.  6 [8,17,29]. 

The  situation  Is  more  complicated  In  the  axisymmetrlc  case  (v*l). 

According  to  Theorem  12,  the  angle  y Is  acute  for  > MQ ( k ) . Therefore, 


A 
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even  In  the  axlsymmetrlc  case  the  first  type  of  domain  of  Influence  (if  it 
exists)  can  hold  only  for  M < M ( k) , when  the  angle  y will  be  oblique, 
which  governs  the  location  of  the  sonic  point  of  the  shock  wave  with  coor- 
dinates R/y,  (R  Is  the  radius  of  curvature  of  the  shock)  relative  to  the 
dividing  curve  ffM^.k),  on  which  y = %tt  [9]  (Fig.  2). 

It  Is  Interesting  to  note  that  domains  of  Influence  of  the  second  type 
(y<it/2)  were  realized  In  computations  of  the  flow  around  a sphere  (k  = 1.4, 
M (k)  = 1.69)  down  to  M =1.1.  The  magnitude  of  the  angle  6 varied  thus  in 
this  case:  6 = u/2  for  = M2  = 3.7  (the  point  K Is  on  the  body);  for 
< M2  the  angle  6 Is  acute  (domain  of  influence  of  the  second  type);  for 
> M2  the  angle  <5  is  oblique  (domain  of  Influence  of  the  third  type). 

Flow  patterns  and  domains  of  influence  for  a sphere  (k=1.4,  and 
M = 1.15,  1.5,  2,3,4),  obtained  from  numerical  computations,  are  presented 

oo 

In  Fig.  7 [17,28]. 

3.3  SUPERSONIC  FLOW  AROUND  BLUNT  BODIES 
WITH  A BREAK  IN  THE  GENERATOR  OF  THE  CONTOUR 

Supersonic  flow  around  blunt  bodies  which  have  a corner  are  calculated 
by  scheme  I or  II  of  the  Dorodnitsyn-Belotserkovskii  method  of  Integral  re- 
lations [8,17,23,24,26]  whose  computational  meshes  are  shown  In  Fig.  8 and 
9. 

The  first  scheme  of  the  Integral  relations  method  is  applied  in  both 
the  domain  up  to  the  corner  and  In  the  supersonic  expansion  around  It,  where 
the  governing  equations  are  written  In  a polar  coordinate  system  centered 
at  the  corner.  The  resulting  system  of  ordinary  differential  equations  is 
Integrated  numerically  along  the  shock  layer  starting  from  the  axis  of  sym- 
metry. In  the  neighborhood  of  the  corner,  where  the  Prandtl  -Meyer  solution 
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holds,  a differential  relation  is  used  which  Is  the  compatibility  condition 
along  the  second  family  of  characteristics  [17]. 


In  the  second  scheme  of  the  method  (Fig.  9),  the  domain  of  Influence  of 
the  bluntness  is  delineated  exactly  because  of  the  simultaneous  construction 
of  the  limiting  characteristic,  which  substantially  Increases  the  accuracy 
of  the  calculation  in  each  approximation.  Here  the  asymptotic  Vaglio- 
-Laurin-Shugaev  solution,  reduced  to  a form  convenient  for  calculations 
[17,25,26],  is  utilized  in  the  neighborhood  of  the  corner  (the  domain  G). 

The  solution  which  describes  plane  and  axisymmetric  transonic  flow  of 
a perfect  gas  in  the  neighborhood  G of  a corner  has  the  form  of  a power 
series  in  the  distance  no  from  the  body  surface  (nQ  = 0)  with  coefficients 
dependent  on  the  corresponding  self-similar  variable  c: 


u ■ 1 + (k+l)"1/3  z u.(c)n  ^2+i^/4,  v = z v.(?)n 
i=0  10  1=0  1 c 

' ■ iwr'Vs/4  ■ % * >v  ■ 

sQ  = (s-s*)/ro*  , 


(3+0/4 


(6) 


where  u ,v  are  velocity  components  (referred  to  the  critical  speed  of  sound) 
along  the  tangent  and  normal  to  the  body  surface  in  the  subsonic  neighbor- 
hood of  the  corner;  quantities  at  the  sonic  point  on  the  body  are  denoted 
with  an  asterisk  (Fig.  8).  The  main  term  of  this  expansion  described  two- 
-dimensional  transonic  potential  flow  in  the  neighborhood  of  the  corner 
(primes  denote  derivatives  with  respect  to  c): 
u0  = g'.  v0=  (79  - 5?g')/4 

where 

g"  - B1(21g  - 25cg')/16,  Bj1  = g'  - (25?2)/16. 

It  is  interesting  to  note  that  the  function  g(c),  corresponding  to  the  poten- 
tial of  this  flow,  can  be  expressed  in  the  parametric  form  [27]: 
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g = (25/42)5_1/8(5Z2+5Z-4)(l-Z)"7/8(l+3Z/5)"9/8C'3  , 

C = -2.5'3/8Z(l-Z)"5/8(l+3Z/5)’3/8C_1 , -5/3 < Z<  1 . ^ 

If  the  scale  factor  C=1 , then  g'(l)=0,  where  g = ( 1 25/56 ) 2-1 /5( -c; ) 7/5  as 
£-►<»  and  g=c3/3  - 675/96  lO-173?^73  as  (refer  to  equation  (8)  below,  for 
the  relationship  between  ? and  Z). 

The  terms  of  the  expansion  (6)  of  a higher  order  of  smallness,  which 
take  into  account  the  rotational  and  axisymmetric  nature  of  the  flow  (u.  , 
v..  ,i  = 1,2,...),  are  found  from  the  solution  of  linear  non-homogeneous 
ordinary  differential  equations.  The  functions  u.j,v..  must  satisfy  the  boun- 
dary conditions  in  both  the  subsonic  domain  (zero  normal  component  of  velo- 
city relative  to  the  body  on  its  surface)  and  in  the  supersonic  domain  where 
the  solution  describes  a flow  of  Prandtl-Meyer  type. 

An  analytic  solution  of  this  system  was  found  in  [17,26]  by  F.V.  Shugaev. 
Let  us  represent  the  quantities  u^,v^  as  the  sum  of  the  particular  solution 

u^J.v/1)  of  the  non-homogeneous  system,  and  the  general  solution 
(2) 

v.j'  ' of  the  homogeneous  system.  If  is  introduced,  so  that  the  solu- 

(2)  (21 

tion  u.'  ,v..v  ' of  the  corresponding  homogeneous  system  would  be  written  as 

d$. 


ui 


(2)  _ d*i 


5 „ “*1 

4 c -ar 


~Si  * vi(2)  = 7~r  *i 
then  a second  order  equation  is  obtained  to  determine 

(16g,-25c2)$."  + [1 6g " + 5(5-2i)c]*i'-(7+i)(3-i)*i  = 0 

Utilizing  the  parametric  representation  (7)  of  the  function  g(c)  and 
making  the  change  of  variables 

ti=*i(l-Z)(7+i)/8(l+  | Z)3(7+i)/8,5  = 1 [2  -*/(3/2(l-Z))]  (8) 

we  obtain  the  hyper geometric  equation 

66(1-5)^"  + (9+2i)(l-2c)t1  - 2(7+1 ) (9+1  )ti  = 0 

Its  solution  is  made  up  of  Jacobi  polynomials 


16 


t,  . [5(H)](5/2+,/3)^)U(l-t)](2i/3+9/2)), 

which  permit  expression  of  the  quantities  u^.v^  in  finite  form. 

Such  a representation  turns  out  to  be  quite  convenient  for  calculations 
in  the  G-domain  since  direct  integration  of  the  original  system  would  give 
rise  to  definite  difficulties.  Computations  have  shown  that  it  is  sufficient 
to  use  the  first  two  or  three  terms  of  the  solution  (6)-(8)  to  construct 
the  flow  in  the  neighborhood  of  an  angular  point;  the  influence  of  the  re- 
maining terms  is  negligible. 

Because  the  limiting  characteristic  of  the  domain  of  influence  of  the 

bluntness  approaches  the  corner  from  the  supersonic  side,  expansion  (6), 

which  is  valid  only  near  the  sonic  line  in  the  supersonic  domain,  turns  out 

to  be  inadequate.  Hence,  it  is  necessary  to  continue  the  solution  Into  the 

purely  supersonic  domain,  which  has  indeed  been  found  as  a correction  in 

ai 

powers  of  r-j  (Fig.  8)  of  the  type  z f.^  J to  the  fundamental  Prandtl-Meyer 
solution. 

The  customary  algorithm  of  scheme  II  is  applied  outside  the  domain  G. 
Systems  of  ordinary  differential  equations  are  integrated  across  the  shock 
layer  between  the  bow  shock  wave  and  the  body,  while  the  additional  condi- 
tions of  "joining"  both  solutions  along  the  boundary  of  this  domain  make 
the  problem  single-valued. 

On  the  whole  we  have  succeeded  In  obtaining  algorithms  which  enable  us 
to  compute  solutions  with  a high  degree  of  accuracy.  The  solution  outside 
. the  zone  of  influence  of  the  bluntness  Is  constructed  by  the  usual  method  of 
characteristics.  Some  results  of  computations  obtained  by  A.  Bulekvaev, 

E.S.  Sedova  and  F.V.  Shugaev  [17,26]  are  presented  in  Figs. 10-12. 

The  geometry  of  the  detached  bow  shock,  sonic  line  and  limiting 
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characteristic  for  perfect  gas  flow  (k=1.4,  Moo=4)  around  a spherical  seg- 
ment with  a sonic  corner  followed  by  a reversed  cone  of  xs30°  half  angle, 
are  shown  in  Fig.  10;  the  corresponding  features  are  shown  in  Fig.  11  for 
x=33°  and  43°30'  with  M^  = 10.  The  variation  of  the  velocity  components 
(along  n,s,  respectively)  on  the  axis  of  symmetry  (i=0),  the  limiting  char- 
acteristic (i=l),  and  the  intermediate  line  (1=2)  Is  given  in  Fig.  12.  It 
is  seen  that  the  behavior  of  u-j  on  the  limiting  characteristic  depends 
strongly  on  the  half-angle  of  the  afterbody. 

3.4  FLOW  AROUND  BLUNT  BODIES  AT  LOW  SUPERSONIC  VELOCITIES 

As  has  already  been  remarked,  the  construction  of  special  numerical 
schemes  would  be  required  for  the  computation  of  mixed  gas  flows  at  low 
supersonic  velocities.  In  this  case,  scheme  III  of  the  method  of  integral 
relations  [17,24,28]  turned  out  to  be  the  most  effective,  wherein  a repre- 
sentation of  functions  In  two  directions  Is  used,  and  the  original  equations 
are  approximated  by  a nonlinear  system  of  algebraic  equations  In  a curvi- 
linear computational  mesh  Fig.  13. 

With  these  schemes  F.D.  Popov  [17]  carried  out  calculations  for  free 
stream  Mach  numbers  down  to  = 1.05  for  the  complete  set  of  gasdynamlcs 
equations  which  Included  effects  of  vortlcity.  The  calculation  of  sonic 
flows  (Mm  = 1,  potential  flow),  where  the  shock  stands  off  at  an  Infinite 
distance,  was  also  carried  out  by  the  method  of  integral  relations  by 
P.I.  Chushkln  [29].  The  results  of  some  of  these  calculations  are  Illustra- 
ted in  Figs.  14-17.  The  geometries  of  the  detached  bow  shock  and  limiting 
characteristic  for  flow  of  a perfect  gas  (k=l .40)  past  a sphere  are  shown 
In  Fig.  14  and  15  for  different  M^.  Experimental  results  of  D.P.  Bedln  and 
G.I.  Mishin  are  indicated  by  the  open  circles.  Fig.  16  shows  the  pressure 
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distribution  along  the  surface  of  a sphere  (referred  to  the  stagnation 
point  pressure)  for  different  M^.  As  the  Mach  number  decreases,  the  pres- 
sure distribution  differs  more  and  more  radically  from  the  Newtonian  dis- 
tribution often  utilized  in  the  approximate  theory  of  supersonic  flow. 

In  Fig.  17d  computational  results  (solid  line)  are  compared  with  the 
experiments  of  V.G.  Maslennikov  et  al . [30]  for  the  bow  shock  standoff  dis- 
tance eQ  along  the  axis  of  symmetry  of  a sphere  as  a function  of  the  free- 
-stream  Mach  number. 

G.M.  Raibinkov  made  a detailed  experimental  investigation  of  the  flow 
behind  a detached  shock  and  presented  the  results  in  the  form  of  tables 
and  graphs  [17].  Figure  17  (a,b,c,e,f)  give  a comparison  between  the 
experimental  data  of  G.M.  Riabinkov  and  the  results  of  calculations  by  the 
method  of  integral  relations  for  a sphere  (6=1)  and  ellipsoids  of  revolution 
(6=0.5  and  2 where  6 is  the  ratio  between  the  vertical  and  horizontal  axes 
of  the  ellipsoid).  The  pressure  distribution  along  the  body,  as  well  as 
the  shape  and  location  of  the  bow  shocks,  are  presented.  It  is  seen  that 
agreement  between  theoretical  and  experimental  results  is  very  good  every- 
where. 
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4.  SECONDARY  SHOCK  WAVES 


4.1  FORMATION  OF  A SHOCK  WAVE  WITHIN  THE  MINIMUM-DOMAIN  OF  INFLUENCE 

It  has  been  shown  in  [6]  that  potential  flow  in  a local  supersonic  zone 
of  definite  type  ("zone  I")  is  spoiled  if  an  arbitrarily  small  length  of  the 
profile  segment  bounding  this  zone  becomes  a straight  line.  Zone  I Is  charac 
terized  by  the  fact  that  both  characteristics  go  from  each  point  of  the  pro- 
file segment  bounding  it  to  the  sonic  line  (see  Fig.  18). 

It  turns  out  that  an  analogous  result  holds  even  in  the  case  of  sym- 
metric rotational  flow  around  a smooth  convex  profile  with  a detached  shock 
[4],  The  existence  of  zone  I here  follows  from  Theorem  13  (Figs.  1,  19). 

Let  us  designate  zone  III  as  a subsonic  domain  of  monotonely  decreasing 
dependence  of  the  entropy  on  the  stream  function.  In  the  case  of  a strictly 
convex  profile,  there  exists  an  e-nelghborhood  of  the  sonic  point  of  the 
profile,  whose  mapping  In  the  hodograph  plane  Is  single-valued. 

THEOREM  19.  A profile  segment  bounding  an  e-nelghborhood  of  the  sonic 
point  of  the  profile  In  zone  III,  whose  mapping  Is  single-valued  in  the  hodo- 
graph plane,  cannot  contain  a straight-line  section  [4], 

Therefore,  making  some  part  of  the  profile  In  the  e-nelghborhood  (for 
X > 1)  Into  a straight  line  results.  Independently  of  its  length,  either  In 
such  a deformation  of  the  sonic  line  or  the  characteristics  so  that  zone  III 
will  not  contain  the  straight  line  section,  or  In  the  formation  of  branch 
lines  which  violate  the  one-sheetedness  of  the  mapping  of  the  orlgnlnal 
e-nelghborhood,  or  In  the  formation  of  a shock  wave.  Apparently,  the  first 
does  not  hold,  since  otherwise  the  solution  of  direct  flow  problem  would  not 
be  unique  or  there  would  not  be  a continuous  dependence  on  the  boundary  con- 
ditions. For  the  direct  flow  problem  It  can  be  shown  that  violation  of  the 


one-sheetedness  of  the  mapping  Is  also  Impossible  because  In  this  case 
changes  would  occur  In  the  flow  which  depend  not  on  the  length  of  the  recti- 
fied section,  but  only  on  Its  location  (under  the  condition  that  the  flow  In 
zone  III  remains  continuous) . These  considerations  lead  to  the  following 
theorem. 

THEOREM  20.  If  the  problem  of  external  flow  around  a profile  with  a 
detached  shock  Is  correct,  then  Its  solution  Is  not  generally  continuous 
(In  the  domain  behind  the  detached  shock  [4]). 


i 


I I 

! 

I I 


4.2  SPOILING  THE  CONTINUOUS  SUPERSONIC  FLOW  DOMAIN  IN  THE  CHARACTERISTIC 
TRIANGLE  BOUNDED  BY  THE  PROFILE  AND  THE  AJOINING  M-DOMAIN 

We  assume  the  existence  of  a smooth  convex  profile  such  that  the  flow  In 
the  region  behind  the  shock  Is  continuous,  and  the  entropy  In  the  M-domain 
monotonely  decreases  with  the  stream  function  (l.e.,  the  shock  has  turned  the 
convexity  towards  the  free  stream).  It  follows  from  Theorem  13  that  the 
boundary  of  the  M-domain  contains  a segment  of  a characteristic  of  the  first 
family  AB,  which  has  one  end  on  the  profile,  and  the  other  on  the  sonic  line 
(Fig.  20). 

Let  A denote  the  triangle  bounded  by  a characteristic  of  the  first 
family  AB,  a characteristic  of  the  second  family  AC,  and  the  profile  segment 
BC.  The  triangle  A adjoins  the  M-domain  from  downstream  (see  Fig.  20). 

THEOREM  21.  Suppose  part  of  a profile  downstream  of  some  point  E Is 
replaced  by  a straight  line  tangent  to  this  point,  and  then  consider  up- 
stream displacement  of  the  point  E.  Before  the  point  E enters  into  the 
M-domain  the  continuous  supersonic  flow  In  the  triangle  A collapses:  either 
a secondary  shock  wave  or  a local  subsonic  zone  arises  In  the  flow  [15]. 
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In  flow  about  a semi -infinite  blunt  wedge  with  a detached  shock  wave, 
either  an  Imbedded  shock  or  a local  subsonic  zone  of  isentropically  retarded 
gas  (Fig.  20)  will  arise  as  the  vertex  angle  Is  gradually  Increased  so  long 
as  the  flow  far  downstream  of  the  bow  shock  remains  supersonic.  An  analogous 
result  is  obtained  in  the  case  of  flow  about  a profile  with  a corner  from 
which  a sonic  line  originates. 

We  shall  restrict  the  proof  of  the  preceding  theorem  to  the  case  of  a 
smooth  convex  profile  in  incident  flows  with  M < M„(k)  (Fig.  20).  In  this 
case  the  first  type  of  minimal  region  of  Influence  of  the  blunt  nose  is  found, 
bounded  In  the  shock  layer  by  the  characteristic  of  the  first  family  AB. 

Integrating  the  compatibility  relations  on  the  characteristics  in  the 
direction  of  increasing  entropy  from  point  A to  the  profile  contour,  we 
obtain  [15] 


eB  - BA  + !1  = BA  + 2lk  * sin2ads  * 

(AB) 

BC  - BA  ‘ *2  = BA  " 2tt  (/B)Sin2adS  ’ 


(9) 


where  a = arc  sin  (1/M),  B is  the  angle  of  inclination  of  the  velocity  vector 
to  the  axis  of  symmetry  and  R Is  the  gas  constant.  Since  both  Integrals  are 
positive  there  must  exist  a point  D on  the  profile  contour,  lying  between 
points  A and  C in  which  Bp  = B^. 

We  subject  the  profile  to  a continuous  deformation,  replacing  part  of 
Its  contour  located  downstream  from  some  point  E,  by  the  tangent  to  the  con- 
tour at  this  point  and  then  displacing  the  point  E upstream  from  its  aft-most 
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position  all  the  way  up  to  the  point  D.  We  shall  show  that  just  as  soon  as 
point  E coincides  with  point  D,  either  one  of  the  aforementioned  assumptions 
is  contradicted  or  the  minimal  region  of  Influence  is  reconstructed.  Indeed, 


1 Inf 
2Rk  E<DF] 


/ s1n2ads 
(AC) 


1 Inf 
3 2Rk  E«[DC] 


f s1n2ads 
(AC) 


f sin2ads 
(AP) 


where  P Is  the  point  of  Intersection  of  characteristics  of  the  second  family 
AC  with  characteristics  of  the  first  family  DP. 

Let  us  now  define  the  point  G on  the  original  profile  by  the  equality 

eG  = eD  " e ' 

In  view  of  the  convexity  of  the  profile,  point  G lies  between  points  D and 
C.  If  point  E Is  located  between  points  D and  G,  l.e..  If 

eD  ^ 6E  > 6G 
then  we  obtain 

ec  ■ be  * 6g  ■ 6d  • c ■ 8a  - ik  eltoc]  ('c)s1n2ads  ’ 

1 .e. , ec  > - I,  which  contradicts  the  second  of  the  inequalities  (9). 

Thus,  while  the  problem  of  determining  supersonic  continuous  flow  In 
the  characteristic  triangle  ABC  from  Cauchy  data  on  the  characteristic  AB 
and  the  Impenetrability  conditions  on  tho  original  profile  BC  had  a solution 
In  the  large,  this  Is  no  longer  true  for  the  profile  deformed  In  the  manner 
Indicated. 

In  an  analogous  way  the  case  M > M (k)  can  be  examined,  as  can  the  case 

oo  0 

In  which  the  sonic  line  emlnates  from  the  corner  on  the  profile.  Results  of 
numerical  computations  agree  very  well  with  the  results  of  the  analytical 
studies  cited  In  this  section  for  smooth  bodies. 
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4.3  FORMATION  OF  IMBEDDED  FLOATING  SHOCK  WAVES  IN  FLOWS  PAST 


PROFILES  WITH  A CORNER 

Let  us  now  examine  [33]  the  flow  about  a profile  with  a convex  corner 
immersed  in  a uniform  supersonic  stream,  especially  when  changes  in  entropy 
on  the  shock  wave  can  be  neglected.  We  shall  also  take  advantage  of  the 
transonic  approximation  to  the  shock  polar  and  the  characteristics  inter- 
secting it. 

First  consider  supersonic  flow  past  a wedge-slab  body  (a  profile  with 
straight-line  segments  OA  and  AF  in  Fig.  21)  with  an  attached  shock  wave. 

In  Fig.  21a,  OB  is  the  straight-line  segment  of  the  shock  wave*,  AB,AD,AC 
are  the  straight-line  segments  of  characteristics  of  the  first  family  while 
BE  is  a characteristic  of  the  second  family.  The  region  behind  the  shock 
wave  in  physical  space  is  mapped  onto  the  velocity  hodograph  plane  n.B  where 
n=(k+l)  (x-1 ) , x is  the  velocity  coefficient, 6 is  the  angle  of  inclination 

of  the  velocity  vector  to  the  axis  of  symmetry,  k is  the  adiabatic  index; 
the  axes  3>y  are  directed  vertically  upward,  the  axes  X,x  horizontally  to  the 
to  the  right.  The  mapping  of  the  region  OBEAO  (Fig.  21b)  is  the  segment  of 
the  characteristic  a^a2  of  the  second  family  e = C - (2/3)n3^;  the  point  a^ 
with  coordinate  B,  lies  on  the  shock  polar  B=n  (n  * \)  (n  -n),  the  point 

I 00  00 

a2  lies  on  the  n axis.  The  equations  of  the  characteristics  and  shock  polar 
are  given  in  the  transonic  approximation. 

Theorem  22.  If  the  flow  field  about  a profile  with  straight  line  seg- 
ments OA  and  AF  behind  an  attached  shock  wave  is  everywhere  supersonic,  then 
It  cannot  be  continuous. 

Let  us  examine  the  fan  of  characteristics  of  the  first  family  emanating 
from  the  corner  point  A (subsequently  denoted  as  fan  A)  In  the  physical  plane 
(Fig.  21a).  If  the  flow  behind  the  shock  wave  is  everywhere  supersonic  and 
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continuous,  then  every  characteristic  of  the  fan  either  intersects  the 
shock  wave  or  extends  out  to  Infinity.  Corresponding  to  this,  in  the  nB 
plane,  each  characteristic  of  the  first  family  passing  through  the  charac- 
teristic a^  must  end  up  either  somewhere  on  the  shock  polar  or  at  the  point 
n,  which  also  lies  on  the  shock  polar  and  represents  the  mapping  of  the 
uniform  rectilinear  flow  at  infinity.  This  represents  a contradiction 
inasmuch  as  characteristics  of  the  first  family  emanating  from  points  of  the 
segment  ca2  in  the  nB  plane  cannot  end  up  on  the  shock  polar;  consequently, 
if  the  flow  is  to  be  supersonic  a secondary  shock  must  arise. 

This  theorem  with  appropriate  changes  in  formulation  can  also  be 
extended  to  the  case  of  a detached  shock  wave. 

It  can  be  shown  [33],  that  the  beginning  of  the  shock  does  not  lie  at 
the  apex  of  the  convex  angle,  but  again  in  the  region  covered  by  the  charac- 
teristics of  the  fan  A.  The  imbedded  floating  shock  is  inclined  downstream 
from  its  origin.  The  end  of  the  shock  of  the  first  family  (where  it  degener- 
ates into  a characteristic)  lies  at  an  infinite  distance  from  the  profile. 

If  a shock  of  the  first  family  intersects  the  main  shock  wave  or  another 
shock,  then  at  the  point  of  intersection  its  intensity  does  not  vanish. 

The  proof  of  these  properties  for  IFS*  is  based  on  the  one-to-one  pro- 
perty of  the  mapping  in  the  hodograph  plane  of  the  apex  of  the  convex  angle 
and  its  neighborhood  covered  by  characteristics  of  the  fan  A. 

Me  denote  by  QR  the  region  bounded  by  the  shock  wave,  the  segment  of 
profile  OA,  the  imbedded  floating  shock,  and  the  last  characteristic  of  the 
fan  A.  In  accordance  with  theorems  of  [33],  the  beginning  of  the  shock  lies 
either  on  this  characteristic  or  on  a characteristic  of  the  first  family 
located  downstream  of  it.  A study  of  the  mapping  of  the  region  QR  in  the 

* IFS  - Imbedded  Floating  Shocks. 
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hodograph  plane  leads  to  the  conclusion  that  the  flow  in  the  region  QR  does 
not  depend  on  the  shape  of  the  profile  over  the  segment  AF.  This  feature 
was  also  observed  in  the  numerical  solutions  of  [17].  In  the  flow  about  a 
profile  with  a straight-line  segment  AF  the  last  characteristic  of  fan  A 
represents  a branch  line. 

Finally,  the  following  theorem  can  be  established: 

THEOREM  23.  In  the  flow  about  a profile  with  a straight-line  segment 
AF  with  = 0 the  imbedded  shock  does  not  intersect  the  bow  shock  wave  at  a 
finite  distance;  apart  from  that,  there  are  no  other  imbedded  floating  shocks 
of  the  first  family. 

In  case  there  exists  a second  shock  in  flow  about  a profile  with  an 
afterbody  segment  AF  having  PQ  > 0,  it  must  intersect  the  bow  shock  wave. 

Up  to  their  point  of  intersection  behind  the  bow  shock  wave,  B < PQ.  The 
bow  shock  wave  after  the  point  of  intersection  with  a second  shock  (or  the 
point  of  intersection  with  the  last  characteristic  of  the  fan  A,  if  there 
is  no  second  shock)  consists  of  an  infinite  number  of  segments  with  curva- 
tures of  opposite  signs*;  the  oscillations  In  the  angle  of  inclination  of 
the  shock  wave  decay  with  Increasing  distance  from  the  profile.  This  result 
also  agrees  well  with  numerical  results  for  bodies  with  corners  [17,26]. 

4.4  CALCULATION  OF  FLOATING  AND  SECONDARY  SHOCK  WAVES 

Imbedded  "floating"  compression  shocks  are  formed  behind  the  bow  shock 
in  the  calculation  of  supersonic  flow  around  blunt  wedges  and  cones  (smooth. 


The  property  of  oscillation  of  the  angle  of  Inclination  of  the  bow 
shock  Is  found  In  connection  with  the  theory  of  propagation  of  dis- 
turbances in  the  case  of  flow  about  wedge-shaped  profiles. 
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with  a break  in  the  generator)  in  definite  flow  regimes  [16,17].  These 
results  demonstrate  that  effects  due  to  viscosity  (i.e.,  boundary 
layer  separation)  do  not  underlie  such  phenomenon  (confirmed  by 
experiment).  In  certain  regimes  of  supersonic  flow  about  blunt  bodies 
(cones,  wedges,  etc.)  imbedded  floating  shocks  (IFS)  arise  in  the  flow  field 
behind  the  bow  shock  which  have  a significant  effect  on  their  aerodynamic 
characteristics.  This  phenomenon  was  also  noted  earlier  in  experiments; 
its  cause,  however,  was  not  particularly  well  understood.  It  was  suggested 
by  M.J.  Lighthill,  for  example,  that  the  shock  was  generated  as  the  result 
of  separation  and  subsequent  reattachment  of  the  boundary  layer  in  the  vi- 
cinity of  the  corner  [33].  The  reasons  for  the  generation  of  IFS,  the 
conditions  under  which  they  arise  (body  configuration,  flight  regime,  etc.), 
and  the  properties  of  such  flows  have  not  been  adequately  studied  to  the 
present  time. 

We  present  here  certain  results  of  computations  related  to  the  genera- 
tion of  IFS.  The  computations  were  carried  out  by  the  aforementioned  me- 
thods within  the  context  of  an  ideal  gas  for  both  smooth  blunt  bodies  and 
profiles  with  a sonic  or  supersonic  corner  [17]. 

Let  us  first  consider  those  properties  which  lead  to  the  generation  of 
IFS.  In  flow  about  blunt  cones  or  wedges  there  occurs  a rapid  deceleration 
of  the  flow  in  passing  from  the  blunt  nose  portion  to  the  straight  portion 
of  the  contour.  This  results  in  crowding  together  of  the  characteristics 
and  in  the  appearance  of  compression  waves  emanating  from  the  body  surface 
behind  the  corner.  In  those  cases  when  the  pressure  gradient  downstream 
from. the  characteristic  bounding  the  region  of  Influence  of  the  blunt  nose 
Is  large  and  the  shock  layer  Is  sufficiently  thick  the  compression  waves 
coalesce  into  a shock  wave  of  initially  zero  Intensity  which  arises  in  the 
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interior  of  the  flow. 

An  imbedded  floating  shock  by  analogy  with  characteristics,  will  be 
Identified  as  a shock  of  the  first  (second)  family  If  the  tangent  to  it  is 
achieved  by  rotating  the  velocity  vector  through  a positive  (negative)  acute 
angle.  Shocks  of  different  families  do  not  connect  smoothly  with  each 
other. 

Every  Instance  of  spontaneous  generation  of  shock  waves  is  always 
associated  with  the  intersection  of  characteristics  of  one  family.  In  such 
flows  there  arises  an  entire  region  doubly  covered  by  characteristics  of  one 
family.  Furthermore,  the  cusp  of  the  envelope  of  these  characteristics 
marks  the  beginning  of  the  shock  wave.  The  cusp  Itself  does  not  yet  belong 
to  the  discontinuity.  The  derivatives  of  velocity,  density,  pressure,  etc., 
become  infinite  at  the  cusp,  and  hence  the  Intensity  of  the  IFS  at  this 
point  must  be  zero  and  a "real"  shock  wave  arises  only  subsequently. 

In  practice  It  is  usually  not  necessary  to  determine  exactly  the  point 
at  which  the  envelope  has  a cusp  or  to  start  out  the  IFS  at  zero  intensity. 

The  shock  wave  Is  considered  Instead  to  emanate  from  the  point  of  Intersec- 
tion of  those  characteristics  of  one  family  which  are  followed  in  the  compu- 
tations. It  should  be  noted  that  IFS,  In  view  of  their  relativity  weak 
Intensity,  can  be  rigorously  constructed  numerically  only  by  direct  applica- 
tion of  the  method  of  characteristics. 

Fig.  22,  taken  from  [34],  shows  examples  of  numerical  results  for  flow 
of  an  Ideal  gas  (k=1.4)  at  = 100  past  a nose  cone  of  semi  vertex  angle 

a>0  = 54.5°  with  conical  afterbodies  of  various  half  angles  u>^  =0°,10°  and  20°. 
The  generator  of  the  body  of  revolution  is  described  by  the  cross-hatched  line 
ACB;  AD  represents  the  bow  shock  wave  and  BD,  the  characteristic  of  the 
second  family  along  which  characteristics  of  the  first  family  Intersect  up 
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to  point  H (If  such  exists)  denoted  by  a dot.  In  Fig.  22c  the  line  CH  is 
the  last  characteristic  of  the  fan  emanating  from  point  C.  The  dot-dash  line 
in  Figs.  22a  and  22b  represents  the  pressure  on  pointed  cones  with  semi 
apex  angles  of  20°  and  10°,  respectively,  for  the  same  Mach  number,  M^  . 
Lengths  are  referred  to  the  radius  of  the  cone  at  the  corner  and  the  surface 
pressures  is  divided  by  the  stagnation  pressure. 

As  indicated  by  the  pressure  distributions,  in  the  third  case,  the  IFS 
is  generated  inside  the  flow  field,  in  the  second,  immediately  behind  the 
bow  shock  wave,  and  in  the  first  case  there  is  no  IFS  in  the  shock  layer. 
Furthermore,  it  should  be  noted  that  as  the  angle  of  the  conical  afterbody 
o>1  is  decreased,  there  is  a corresponding  increase  in  the  pressure  gradient  at 
the  corner,  and  the  shock  layer  becomes  relatively  thicker. 

The  calculations  were  obtained  by  the  method  of  characteristics  for  an 
ideal  gas;  the  initial  conditions  for  the  calculations  were  obtained  from 
the  exact  solution  for  a pointed  cone,  insuring  a high  degree  of  accuracy 
in  the  computations.  These  examples  show  that  the  generation  of  IFS  is  not 
dependent  on  flow  properties  related  to  viscosity  (i.e.,  separation  of  the 
boundary  layer)  and  moreover  that  their  occurrance  is  not  related  to  the  use 
of  "coarse"  initial  data  for  the  numerical  calculation  of  the  supersonic 
region  by  the  method  of  characteristics. 

Thus  the  regimes  of  flow  about  blunt  cones  and  wedges  which  are  condu- 
cive to  the  development  of  a positive  pressure  gradient  behind  the  corner 
and  a sufficiently  thick  shock  layer  lead  to  the  generation  of  IFS,  the  posi- 
tion of  which  is  determined  by  the  intersection  of  two  neighboring  charac- 
teristics of  the  same  (first)  family. 

We  shall  now  consider  several  more  examples.  Shock  waves  (Fig.  23a) 
and  surface  pressure  distributions  (Fig.  23b)  are  shown  for  plane  (v=0)  and 
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axisymmetric  (v=l)  bodies  described  by  a circular  nose  with  a corner  at 
x=68°  and  w=0,  while  Fig.  24  shows  the  pressure  distribution  along  smooth 
circular  cones  (solid  line)  and  wedges  (dotted  line)  for  various  apex 
angles  [8,16].  From  these  figures  it  is  apparent  that,  other  things  being 
equal  (the  flow  regime,  contour  shape,  etc.),  the  plane  blunt  body,  in  con- 
trast to  the  axisymmetric  (or  a body  with  a sonic  corner  in  contrast  to  a 
smooth  one)  Introduces  a greater  disturbance  into  supersonic  flow,  and  the 
resulting  shock  layer  is  thicker.  Furthermore,  in  the  case  of  blunt  plates 
u=0  and  wedges  to^O  (plane  flows),  and  also  for  bodies  with  a sonic  corner 
(Fig.  25)  a zone  of  overexpansion  arises  immediately  after  the  blunt  nose 
followed  at  the  shoulder  by  a positive  pressure  gradient  which  Increases 
rapidly  with  increasing  x.  It  follows  that,  all  other  things  being  equal, 
conditions  more  favorable  to  the  generation  of  IFS  are  realized  in  the  plane 
two-dimensional  case  or  for  bodies  with  a sonic  corner. 

Fig.  26  shows  the  main  bow  shock  wave  and  IFS  which  arise  in  the  super- 
sonic zone  in  the  case  of  flow  about  spherically  blunted  cones,  which  have 
a sonic  corner  Cx=30°;  k=l .4)  [17].  Results  are  shown  for  cones  with  semi 
vertex  angles  of  u)=-5°,  0°,  10°  and  Incident  Mach  numbers  of  M =4  (the 
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solid  curve)  and  M^  = 6 (dotted)  curve.  The  intensity  of  the  imbedded 
floating  shock  first  rapidly  increases  and  then  gradually  decreases  with 
Increasing  distance  from  the  blunt  nose.  In  the  cases  cited  above  the  angle 
through  which  the  flow  passing  through  the  imbedded  shock  is  deflected 
achieves  a maximum  of  only  a few  degrees,  while  at  a distance  of  some  30  to 
40  nose  radii  this  angle  amounted  to  only  a few  seconds.  V.F.  Ivanov  first 
computed  tables  of  the  locations  of  floating  shocks  [17].  It  Is  interesting 
to  note  that  when  floating  shocks  formed,  the  shape  and  location  of  the  main 
shock  waves  (in  the  region  lying  ahead  of  their  intersection  with  secondary 
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shocks)  coincided  (for  identical  values  of  incident  Mach  number)  for  various 
half-angles  of  the  cone  in  the  range  -10°  < u < +10°  (Fig.  27).  This  pro- 
perty is  explained,  apparently,  by  the  substantial  extent  of  the  region  of 
influence  of  a blunt  nose,  which,  for  such  regimes  of  flow,  includes  a rela- 
tively large  portion  of  the  main  shock  wave.  Fig.  28  illustrates  the  effect 
of  the  corner  point  on  the  position  of  the  limiting  characteristics. 

An  analogous  situation  regarding  the  region  of  influence  is  also  ob- 
served for  smooth  blunt  cones.  Fig.  29  shows  the  coordinates  of  the  shock 
wave  for  smooth  cones  with  spherical  noses.  The  coordinates  of  the  point  B, 
through  which  passes  the  characteristic  bounding  the  region  of  influence 
of  the  nose  are  shown  for  various  cone  half-angles  u>. 

Fig.  30  shows  examples  of  flow  with  IFS  obtained  by  V.I.  Kosarev  for 
axisymmetric  double-cone  bodies  with  a = 0°  and  5°  at  M^  = 6 ( k=l . 4) ; these 
conditions  insure  isentropic  compression  of  the  flow  at  the  juncture  between 
the  cones.  Flow  patterns  (IFS  are  shown  by  dotted  lines)  and  longitudinal 
surface  pressure  distributions  are  shown  for  the  planes  <j>  = 0,  ir/2,  it.  Here, 
too,  the  generation  of  IFS  is  clearly  seen  when  the  conditions  of  a suffi- 
ciently thick  shock  layer  and  positive  pressure  gradient  are  present  (a  = 0°; 
a = 5°  at  4>  = 0 , ir/2). 
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5.  LOCAL  SUPERSONIC  ZONES 


Recently,  systematic  calculations  of  "supercritical"  and  transonic  ro- 
tational flows  have  been  carried  out  by  the  Computing  Center  of  the  Academy 
of  Sciences  of  the  U.S.S.R.  [37,38].  These  apply  to  plane  and  axisymmetrlc 
bodies  and  are  based  on  the  method  of  "large  particles"  [31,32,37].  In  this 
paper,  without  discussing  the  details  of  the  method,  we  shall  only  be  con- 
cerned with  the  formulation  of  the  boundary  conditions  of  the  problem  for 
this  class  of  flows  and  present  some  computational  results. 

The  development  of  the  unsteady  "large  particles"  method  was  stimulated 
by  the  works  of  F.  Harlow  [35]  and  others.  We  adhere  to  the  "particle-in- 
-cell"  method.  Nevertheless,  it  seems  reasonable  that,  for  gas  dynamics 
problems,  we  should  not  confine  ourselves  to  the  discrete  model  of  a con- 
tinuous medium  comprising  a combination  of  particles  of  a cell  of  fixed  mass. 
Instead,  we  consider  continuous  flows  of  "large  particles"  whose  mass  coin- 
cides with  the  mass  of  an  Ruler  cell  at  a given  instant  of  time.  In  this 
formulation  each  time  cycle  is  divided  into  three  stages  ("Eulerian",  "La- 
grangian",  and  "Final"),  with  greater  attention  given  to  the  development  of 
a numerical  algorithm  which  might  be  used  for  a wide  class  of  problems  con- 
cerned with  motion  of  a compressible  gas  [31,32,35-38]. 

As  a result,  we  have  obtained  some  divergence  form  (conservative),  dis- 
sipative-steady difference  schemes  which  allow  us  to  consider  a wide  class 
of  transonic  problems  as  well.  Developments  In  this  direction  have  been 
carried  out  since  1965. 

5,1  "LARGE  PARTICLES"  METHOD 

Let  us  briefly  describe  the  main  principles  of  the  "large  particles" 
method.  The  region  of  integration  is  covered  by  a fixed  (over  space)  Euler 
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net  composed  of  rectangular  cells  with  sides  Ax,  Ay  (or  az,  Ar  in  a cylindri- 
cal coordinate  system). 

In  the  first  ("Eulerian")  stage  of  calculations  only  those  quantities 
change  which  are  related  to  a cell  as  a whole,  and  the  fluid  is  supposed  to 
be  momentarily  decelerated.  Hence,  the  convective  terms  of  the  form 
div(<t>pw),  where  $ = (1  ,u ,v,E) , corresponding  to  displacement  effects,  are  omit- 
ted in  the  governing  (initial)  equations.  Then  it  follows  from  the  equation 
of  continuity,  in  particular,  that  the  density  field  will  be  "frozen"  and  the 
initial  system  of  equations  will  be  of  the  form 

(5.1) 


plr  + lf=0*  p!t+!£=0’  p ft  + div(ew)  = 0 


Here  we  have  used  both  the  simplest  finite-difference  approximations 
and,  to  improve  the  calculation  stability,  the  schemes  of  the  method  of  inte- 
gral relations  [36],  in  which  "sweeping-through"  approximations  of  the  inte- 
grands with  respect  to  rays  (N  = 3,4,5)  are  used  [37]. 

In  the  second  ("Lagrangian" ) stage  we  find  mass  flows  AMn  across  the 
cell  boundaries  at  time  tn  + At.  At  this  stage  we  assume  the  total  mass  to 
be  transferred  only  by  a velocity  component  normal  to  the  boundary.  Thus, 
for  instance, 

AMi+*s, j = (pi+5s,j)  (ui+is,j)  Ay  At  ’ etc’ 

The  brackets  ( ^ denote  the  values  of  p and  u across  the  cell  boundary.  The 
choice  of  these  values  is  extremely  important  since  they  substantially  influ- 
ence the  stability  and  accuracy  of  the  calculations.  Consideration  of  the 
flow  direction  is  essential  to  the  various  possible  ways  of  writing  down 
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First  and  second  order  accurate  representations  of  AMn  are  considered. 
These  are  based  on  central  differences,  without  account  being  taken  of  the 
flow  direction,  as  well  as  by  means  of  the  discrete  model  of  a continuous 
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i medium  comprising  a combination  of  particles  of  a fixed  mass  in  a cell. 

Lastly,  in  the  third  ("Final")  stage  we  estimate  the  final  fields  of 
the  Euler  flow  parameters  at  the  instant  of  time  tn+1  = tn  + At  (all  the  er- 

r 

l rors  in  the  solution  are  "removed").  As  was  pointed  out,  the  equations  at 

■f 

this  stage  are  laws  of  conservation  of  mass  M,  momentum  P and  total  energy 
E written  down  for  a particular  cell  in  the  difference  form 

Fn+1  = Fn  + EAFnbdry  where  F * (M,P,E)  (5.3) 

According  to  these  equations.  Inside  the  flow  field  there  are  no  sources 

■f 

or  sinks  of  M,  P and  E and  their  variations  In  time  At  is  caused  only  by  in- 
teraction at  the  external  boundary  of  the  flow  region. 

It  follows  from  the  very  character  of  the  construction  of  the  calcula- 
tion scheme  that  a complete  system  of  nonstationary  gas  dynamics  equations  is 
essentially  solved  here,  while  each  calculation  cycle  represents  a completed 
process  in  a given  time  Interval.  Consequently,  the  governing  nonstationary 
equations,  subject  to  the  boundary  conditions  of  the  problem,  are  satisfied 
so  that  the  real  fluid  flow  at  the  time  in  question  is  determined. 

Thus,  the  "large  particles"  method  allows  us  to  obtain  the  characteris- 
tics of  nonstationary  gas  flows  and  as  a consequence  of  their  stability,  the 
asymptotic  steady  state  as  well.  Such  an  approach  is  especially  applicable 
to  problems  in  which  a complete  or  partial  development  of  physical  phenomena 
with  respect  to  time  takes  place.  For  example.  In  studying  transonic  gas 
flows  and  flows  around  finite  bodies,  flow  In  local  supersonic  zones  and  se- 
paration regions  develop  comparatively  slowly  while  the  major  part  of  the 
field  develops  rather  rapidly.  Our  Investigation  Is  wholly  devoted  to  sys- 
tematic calculations  of  a wide  class  of  compressible  flows  Involving  tran- 
sonic regime,  discontinuities,  separation  and  "injection". 


The  divergence  forms  of  the  differential  and  difference  equations  are 


considered  in  the  "large  particles"  method;  the  energy  relation  for  total 
energy  E is  used;  different  kinds  of  approximations  are  used  in  the  1st  and 
2nd  stages;  additional  density  calculations  are  introduced  in  the  final 
stage,  which  helps  to  remove  fluctuations  and  makes  it  possible  to  obtain 
satisfactory  results  with  a relatively  small  network  (usually  1-2.5  thousand 
cells  are  used).  All  this  results  In  completely  conservative  schemes,  i.e., 
laws  of  conservation  for  the  whole  net  region  are  an  algebraic  consequence  of 
the  difference  equations.  Fractional  cells  are  introduced  for  the  calcula- 
tion of  bodies  with  a curvature  in  the  slope  of  the  contour  [37]. 

The  investigation  of  these  schemes  (approximation  problems,  viscosity, 
stability,  etc.)  was  carried  out  successively  for  the  zero,  the  first  and 
the  second  differential  approximations  [31,32,37].  These  investigations 
show  that  the  "large  particles"  method  yields  divergence-conservative  and 
dissipative-steady  schemes  for  "sweeping- through"  calculations.  This  enables 
us  to  carry  out  stable  calculations  for  a wide  class  of  gas  dynamics  problems 
without  Introducing  explicit  terms  with  artificial  viscosity.  It  may  be  of 
particular  significance  In  studying  flows  around  bodies  with  a curvature  in 
the  slope  of  the  contour  since  the  ways  of  introducing  explicit  terms  with 
artificial  viscosity  are  different  for  whole  and  fractional  cells.  Moreover, 
by  varying  only  the  second  stage  of  the  calculation  procedure  we  can  arrive 
at  the  conservative  "partlcle-ln-cell"  method  so  that  the  calculational 
algorithm  Is  of  general  use. 

As  for  discontinuities  the  approximate  viscosity  In  the  scheme  (dissi- 
pative terms  In  difference  equations)  results  in  stable  calculations  with  a 
"smearing"  of  shock  waves  over  several  computational  cells  and  the  formation 
of  a thick  boundary  layer  near  the  body.  It  should  be  stressed  that  the 
magnitude  of  the  approximate  viscosity  is  proportional  to  a local  flow 
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velocity  and  to  the  dimensions  of  the  difference  net,  therefore  its  effect 
is  practically  evident  only  in  zones  with  large  gradients. 

The  boundary  conditions  of  the  problem  were  realized  by  introducing 
layers  of  "fictitious"  cells  along  the  region  boundaries  [31,32,37].  Condi- 
tions on  the  body  closely  approximate  the  "no-slip"  condition.  For  example, 
in  Fig.  31  where  the  density  profile  is  given  for  the  conditions  of  "no-slip" 
(dashed  line)  and  vanishing  of  the  normal  component  of  velocity,  identified 
"non-flow"  (solid  line),  even  in  the  vicinity  of  the  body  the  difference 
between  these  two  cases  is  insignificant  and  at  some  distance  from  the  body 
it  disappears  entirely. 

It  turned  out  that  the  right  hand  "open"  boundary  of  the  region  intro- 
duces the  greatest  distrubances  in  the  calculation  of  transonic  gas  flows. 

To  evaluate  its  influence  and  determine  the  optimum  dimensions  of  the  net, 
the  calculations  were  carried  out  for  nets  of  different  sizes;  "matching"  of 
the  flow  fields  was  enforced  by  utilizing  one  of  the  internal  columns  as  the 
initial  one  for  a new  field.  Finally  a comparison  between  the  asymptotic 
(steady)  state  and  experiment  was  made  as  well  [31,32,37,38]. 

Fig.  32  shows  the  results  of  the  calculations  of  a "super-critical"  flow 
(M^  = 0.9)  over  regions  of  various  extent,  t/R  - 2 to  3 in  the  flow  around  a 
semi-finite  cylinder  (as  if  the  body  were  "moved  into"  the  stream).  If  the 
flow  field  ahead  of  the  body  is  established  rather  quickly,  then  the  flow 
downstream  from  the  corner  become  steady  only  past  1/ R * 2 to  3.  The  results 
of  the  calculations  with  a coarse  net  are  shown  in  Fig.  32d  while  the  region 
of  Figs.  32a-32c  (which  utilize  about  2.5  thousand  nodes)  is  indicated  by  a 
dashed  line. 

A comparison  Is  given  in  Fig.  33  between  the  results  of  the  calculations 
(solid  lines)  and  analytical  data  for  the  asymptotic  steady  state  (dashed 
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line)  for  a sonic  flow  around  the  same  body  [31,32].  In  this  figure  "1"  is 
the  sonic  line,  "2"  the  limiting  characteristic  curve,  "3"  is  the  line  show- 
ing the  departure  of  the  velocity  vector  from  the  horizontal,  and  "4"  is  the 
shock  wave.  It  should  be  noted  that  good  agreement  is  observed  with  the 
analytic  data  already  at  a distance  of  2-3  radii  from  the  body. 

5.2  CALCULATION  OF  TRANSONIC  AND  "SUPERCRITICAL"  FLOWS 

Some  computational  results  obtained  by  the  "large  particles"  method 
[37,38]  for  transonic  flows  around  two-dimensional  and  axisymmetric  bodies 
will  now  be  described.  For  purposes  of  this  discussion  the  supercritical 
regimes  of  transonic  flows  around  bodies  will  be  characterized  by  the  value 
of  the  critical  Mach  number  of  the  oncoming  flow  M*  (i.e.,  when  a sonic 
point  first  develops  on  the  body)  as  well  as  by  the  extent  of  the  local 
supersonic  zone  (as  compared  to  a characteristic  dimension  of  the  body)  and 
by  its  intensity  (say  the  maximum  supersonic  velocity  M realized  in  the 
zone) . 

Fig.  34  (series  1 - 8)  presents  the  flow  field  patterns  (lines 
M = constant)  for  a 24 % circular  arc  profile  (v  = 0)  extending  from  purely 
subsonic  (M^  = 0.6)  to  supersonic  regimes  (Mro  = 1.5).  Successive  flow  fields 
for  increasing  M^  depict  transition  through  the  critical  Mach  number  (here 
M*  = 0.65)  and  the  formation  and  development  of  a local  supersonic  zone. 

The  supercritical  flow  around  this  profile  is  observed  for  0.7  < Ma><l 
(Figs.  34  :2-7).  One  can  distinctly  see  the  position  of  the  shock  In  the 
region  of  crowded  M = constant  lines  which,  together  with  the  sonic  line, 
bound  the  local  supersonic  flow.  The  region  of  subsonic  velocities  is  located 
behind  the  shock  wave.  When  the  velocity  of  the  oncoming  flow  Increases, 
the  flow  disturbances  produced  by  the  body  die  out  at  a large  distance  from 
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the  body.  With  M > 0.9  the  zone  becomes  considerable  both  in  size  and  in 
intensity  (supersonic  velocities  up  to  M = 1 .7  to  1.8  are  attained)  and  in 
case  of  a sonic  flow  (Fig.  34.7)  lines  of  the  level  M = 1 extend  to  infinity. 

The  asymmetry  of  the  whole  flow  pattern  is  noticeable  (even  at  purely 
subsonic  velocities  - Fig.  34.1)  which  results  from  non-potentiality  of  the 
flow  (supercritical  regimes)  and  from  the  presence  of  viscous  effects  as 
well  (formation  of  a wake  behind  the  body). 

In  the  case  of  supersonic  flow  past  this  profile  (Fig.  34.8,  = 1.5) 

a shock  wave  developes  ahead  of  the  body  which  bounds  the  disturbed  region. 
Behind  the  wave,  subsonic  velocities  occur  in  the  vicinity  of  the  axis  of 
symmetry  away  from  which  the  flow  velocity  along  the  contour  of  the  body  in- 
creases and,  as  a result,  a "terminal"  shock  occurs  near  the  stern  of  the 
body. 

For  comparison  the  results  of  calculations  by  the  above  method  for  flow 
around  a 24%  axi symmetric  "spindlelike"  body  (v  = 1)  are  given  in  Fig.  35 
with  0.8  < Mro  < 2.5.  In  this  case  a critical  regime  already  occurs  at 
M*  = 0.86;  local  supersonic  zones  as  compared  to  the  plane  case  are  less  de- 
veloped and  of  weaker  intensity  (for  example,  values  of  M * 1.3-1 .4  are 
realized),  although,  naturally,  the  main  features  of  a transonic  flow  are 
quite  evident  (see  also  [38]). 

In  Fig.  36  a comparison  is  given  between  the  flow  fields  calculated  by 
the  above  method  (solid  line)  and  those  of  the  Wood  and  Gooderum  experiment 
(dashed  line)  for  subcrltlcal  (Fig.  36a,  M = 0.725)  and  supercritical  (Fig. 
.36b,  M*  = 0.761)  flows  around  a 12%  profile  (results  of  both  the  calculations 
and  the  experiment  Indicate  M*  = 0.74). 

As  compared  to  the  corresponding  plane  cases  with  the  same  value  of  6 
It  can  be  seen  that,  in  axisymmetrlc  flows,  an  increase  Is  seen  in  the  critical 
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Mach  number  (with  6 = 1 2%,  from  0.74  to  0.89,  and  from  0.65  to  0.86  for  6 = 24%). 

In  axisymmetric  flow  local  supersonic  zones  are  of  smaller  size  and  lower 
intensity,  for  instance,  the  maximum  Mach  number  in  these  zones  is  M = 1.2 
for  v = 1 as  compared  to  M = 1 .5  for  v = 0 when  6 = 12%,  and  M = 1.3:  1.7  when 
6 = 24%.  However  greater  asymmetry  of  the  flow  picture  about  the  vertical 
axis  is  observed  for  axisymmetric  bodies;  behind  the  body  we  find  a stronger 
wake,  which  exerts  quite  a marked  influence  on  the  flow  picture,  and  finally, 
the  shock  layer  is  considerably  narrower  in  supersonic  axisymmetric  flow. 

Analysis  of  internal  check  tests  as  well  as  the  results  of  various  com- 
parisons indicate  that  the  computational  error  of  the  "large  particles"  me- 
thod usually  does  not  exceed  several  per  cent.  These  calculations  were 
carried  out  using  a Soviet,  BESM-6  computer;  the  time  of  the  calculation  in 
this  case  did  not  exceed  an  hour. 

Some  numerical  results  obtained  by  the  large  particle  method  for  flow 
past  more  complicated  types  of  bodies  will  now  be  considered.  In  Fig.  32  we 
have  plotted  the  M = constant  lines  for  supercritical  flow  past  a semi-infi- 
nite cylindrical  block  (here,  M*  = 0.70,  which  is  in  good  agreement  with  the 
experimental  value  M*  = 0.69  obtained  by  Stanbrook  [40])  while  Fig.  37  refers  to 
the  case  of  sonic  flow  past  a space  vehicle  of  the  Apollo  type.  The  stream- 
lines, bow  shock  wave  and  sonic  line  in  supersonic  flow  past  a short  cylin- 
drical body  and  a sphere  are  shown  in  Fig.  38.  These  numerical  results 
clearly  indicate  the  formation,  behind  "badly  streamlined"  bodies  (Figs.  37 
and  38),  of  closed  recirculation  zones,  localized  behind  the  tail,  which  are 
isolated  from  the  external  flow  by  a contact  surface  (marked  by  a dashed 
line  in  Fig.  38).  The  flow  in  the  recirculation  zones  is  strongly  subsonic 
and  extremely  rarefied  (i.e.,  both  the  gas  density  and  pressure  are  small  In 
them) . 
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An  important  pratical  application  of  the  large  partical  method  is  de- 
scribed by  near  sonic  flow  (M^  = 1)  past  a body  that  has  a jet  (with  para- 
meters M , p , u , v ) "blown  out"  opposite  to  the  main  flow,  thus  changing 
the  entire  flow  field.  In  this  case,  auxiliary  conditions  for  the  jet 
flow  must  be  prescribed,  usually  sonic  velocity,  with  the  jet  parameters 
specified  in  appropriate  cells  on  the  body  surface.  The  field  of  constant 
Mach  number  contours  are  shown  in  Fig.  39a,  for  sonic  flow  (M^  = 1.0)  past 
a sphere  which  has  an  axial  sonic  jet  (M  = 1.0;  p = 2.9,  u = 1.0,  v = 0) , 

v t U V 

issuing  upstream  from  a nozzle  located  on  its  axis  of  symmetry;  Fig.  39b 
gives  corresponding  results  for  supercritical  flow  at  = 0.9.  In  the  pre- 
sence of  a jet  the  flow  around  the  body  becomes  much  more  complicated  and  is 
unsteady  in  the  mixing  region.  The  flow  field  splits  into  two  sub-regions; 
the  gas  coming  from  the  jet  passes  through  one  of  them,  and  the  external  flow 
through  the  other.  The  trend  towards  this  division  may  be  seen  in  Fig.  39a, 
and  becomes  even  more  evident  in  the  supercritical  case  (Fig.  39b).  This  de- 
marcation is  given  by  the  contact  surface  (shown  by  the  dashed  line  in  Fig. 
39b)  which  separates  the  external  flow  from  the  injected  gas.  A separated 
wake  is  also  formed  behind  the  body. 
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6.  ON  THE  QUESTION  OF  EXTREMALITY  OF  THE  ENTROPY  ON 


THE  STAGNATION  STREAMLINE  IN  THREE-DIMENSIONAL  FLOW 

Of  great  interest  in  three-dimensional  supersonic  flow  is  the  question 
as  to  whether  the  stagnation  streamline  (streamline  arriving  at  the  body) 
passes  through  a point  on  the  shock  where  it  is  orthogonal  to  the  velocity 
vector. 

It  has  been  shown  in  [18]  that  if  the  entropy  does  not  have  a maximum 
on  the  surface  of  an  axisymmetric  body  in  supersonic  three-dimensional  flow, 
then  the  vorticity  on  the  body  surface  is  a discontinuous  function  of  the 
angle  of  attack.  This  property  was  found  under  the  assumption  that  the  cur- 
vature of  the  shock  was  finite  at  its  point  of  intersection  with  the  axis  of 
symmetry  (it  was  later  established  that  this  assumption  is  always  true). 

An  investigation  of  the  local  three-dimensional  flow  in  the  neighborhood 
of  the  stagnation  point  on  a body  showed  that  the  vorticity  is  generally  not 
bounded  on  the  body  surface  [19].  In  the  general  case  of  flow  around  a body 
with  a detached  shock,  an  analogous  result  has  been  obtained  under  the  as- 
sumption of  the  existence  of  a stagnation  streamline  in  [20]. 

THEOREM  24.  For  the  normal  derivative  of  the  velocity  on  a body  surface 

to  be  bounded,  it  is  necessary  that 
conditions  for  the  existence  of  an  ex- 
tremum in  the  entropy  be  satisfied  on 
the  stagnation  streamline. 


Figure  40.  Special  coordinate  system 
generated  by  the  stagnation 
streamline. 


If  these  conditions  are  not  satisfied,  then  the  body  surface  is  almost 
everywhere  the  envelope  of  surfaces  of  equal  velocity  magnitude;  a weak  dis- 
continuity is  propagated  along  two  body  streamlines  which  emerge  from  the 
stagnation  point  (Fig.  40).  The  construction  of  an  appropriate  example  in 
[20]  showed  that  the  extremal ity  of  the  entropy  on  the  stagnation  streamline 
is  generally  insufficient  for  the  boundedness  of  the  normal  velocity  derivative 
on  the  body  surface. 

It  should  also  be  noted  that  an  approximate  solution  of  the  problem  of 
three-dimensional  hypersonic  flow  around  a body  [21]  has  shown  that  the  en- 
tropy does  not  take  on  a maximum  value  on  the  body  surface.  An  analogous 
deduction  has  been  made  in  [22]  from  an  experimental  investigation  of  the 
problem.  Numerous  calculations  ([17]  et  al.)  of  the  flow  around  axisymme- 
tric  bodies  at  angles  of  attack  over  a wide  range  of  flow  conditions 
(1  < < 20,  0 < a < 25°)  have  shown  that  the  difference  in  entropy  on  the 

body  surface  from  the  maximum  value  is  of  the  same  order  as  the  computational 
error. 

Therefore,  although  the  connection  between  the  extremum  of  the  entropy 
on  the  stagnation  streamline  and  the  flow  properties  at  the  body  surface  has 
indeed  been  established  in  the  papers  cited,  the  question  posed  herein  re- 
mains open  for  a rigorous  theoretical  analysis. 
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Fig.  7.  Domains  of  influence  in  flows  of  a perfect  gas  (k  * 1.4)  around  a 
sphere:  II  for  1.1  < < M2  = 3.7;  III  for  3.15  < < « . 
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Fig.  8.  Computational  mesh  for  Scheme  1 of  Fig.  9.  Computational  mesh  for 

the  Method  of  Integral  Relations.  Scheme  II  of  the  Method 

of  Integral  Relations. 


Fig.  10.  Geometry  of  supersonic  flow 
with  M^-4  past  a spheri- 
cally-capped reversed  cone 
of  x * 30°  half  angle. 


Fig.  11.  Geometry  of  flow  around 

spherically-capped  noses  of 
reversed  cones  with  M » 10 

and  x - 33°  and  43°30’ . 


Fig.  12.  Distribution  of  the  normal,  u^,  and  tangential,  v^  velocity 

components  across  the  shock  layer  of  the  bodies  of  Fig.  11. 
calculated  from  Scheme  II,  along  the  axis  of  symnetry  i-0, 
along  the  limiting  characteristic,  i » 1 and  the  intermediate 
line,  i » 2. 
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Fig.  13.  Computational  me sh  for 
Scheme  III  of  the  Method 
Integral  Relations. 


Fig.  14.  Geometry  of  Detached  Bow  Shock  in  super 
sonic  flow  of  a perfect  gas  (k  ■ 1.4) 
past  a sphere. 


P.(i) 


Fig.  IS.  Comparison  of  nunerical 
results  with  experiment 
(points)  of  Bedin  and 
Mishin  [30]  for  super- 
sonic flow  (k«1.4)  past 
a sphere. 


Fig.  16.  Calculated  surface  pressure  dis 
tribution  along  a sphere  with 
k - 1.4. 
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Fig.  18.  Local  supersonic  "Zone  I"  in  a 
subsonic  potential  flow  accord- 
ing to  Nikol'skii  and  Taganov 


(a) 


^*8*  Ca)  Zone  III  in  the  physical  plar 
Cb)  Mapping  of  Zone  III  in  the  he 
graph  plane. 


Fig.  20.  The  collapse  of  the  continuous  supersonic  flow  in  the  characteristic 
triangle  A which  adjoins  the  minimal  domain  of  influence  as  the  angle 
of  the  infinite  blunted  wedge  increases. 


Fig.  21.  Supersonic  flow  past  a wedge-slab  body  with  attached  shock  wave:  AB,  AC,  AD 
characteristics  of  the  first  family,  BE  of  the  second  family. 
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circular  nose  and  comer  at  x = 68°. 


Bow  shock  and  imbedded  floating 
shock  for  spherically  blunted 
bodies  with  a sonic  comer. 
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Fig.  30  Imbedded  floating  shock  and  sur- 
face pressure  distribution  over 
a double  cone  in  supersonic  flow, 

M * 6,  k=  1.4  with  a » 0°  and  5°. 


Fig.  31  Density  profiles  for  boundary 
conditions  of  no-slip  ( — ) 
and  vanishing  of  velocity  com 
ponent  normal  to  surface ( ) 


Fig.  32  Lines  of  constant  Mach  nunber  (isotachs)  in 

supercritical  flow  (M^'0.9)  around  flat-faced 

cylinders  of  various  lengths:  (a)  Z/R«0.56, 
Cb)  l/R-2.0,  (c)  t/R-2.72,  (d)  t/R-7. 


Fig.  33  Comparison  of  nunerical  results 
( — ) with  the  asymptotic  steady 
state  [31 ,32]  for  sonic  flow 
past  a flat-faced  cylinder:  1 - sonic  line 
2 - limiting  characteristics,  3-  locus  of 
horizontal  velocity  vector,  4 - shock  wave 


Fig.  35  Isotachs  in  transonic  flow  around  a 24  per  cent  axisymmetric  body 

(generated  by  revolution  of  a circular  arc  profile) ; critical  Mach 

number,  M * = 0.86. 

7 00 


Fig.  36  Isotachs  in  transonic  flow  around  a 12  per  cent  circular  arc  airfoil;  calculated  by 
the  "large  particles"  method , experiment  of  Wood  and  Go ode rum  [39] . 
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This  report  considers  some  problems  of  transonic  gasdynamics  related 
to  the  theory  of  planar  and  three-dimensional  flows  of  a perfect  gas  which 
have  been  worked  out  mainly  by  the  author  and  his  collaborators.  Results 
of  analytical  investigations  of  a number  of  exact  properties  of  solutions 
which  describe  the  flow  around  bodies  with  p rfptarhpd  <;horkp  obtained  by 
E.G.  Shifrin,  O.M.  Belotserkovskii,  et  al.?are  presented  without  proof  in 
the  form  of  theorems  (proofs  are  to  be  found  in  the  papers  cited).  Results. 


DD 


EDITION  OF  I NOV  IS  1$  OBSOLETE 

S/N  0102-  LF-0I4-460I 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whwi  Omit  I 


SECURITY  CLASSIFICATION  OF  This  RACE  (Whmn  Dim  Bntmrud) 
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^of  numerical  solutionsi (due  to  O.M.  Belotserkovskii,  Yu  M.  Davidov,  F.D.  Popov, 
et  al . )*?  obtained  by  the  methods  of  integral  relations  and  £T.arge  particles  , 
are  presented  for  the  complete  problem  of  transonic  flow  around  blunt  bodies 
whose  profiles  possess  sharp  comers.  The  boundaries  of  distinct  minimal 
domains  of  influence  in  mixed  flow  about  blunt  bodies  are  also  examined.  __ 
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